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Q ■ Abstract. We are going to review recent advances in the theory of exclusive nonlep- 

D ' tonic B decays. The emphasis is going to be on the factorization hypothesis and the role 

^0 , of nonfactorizable contributions for nonleptonic B decays. In particular, we will discuss 

■ more in detail calculations of nonfactorizable contributions in the QCD light-cone sum 

' rule approach and their implications to the B —* nn and B — > J/ipK decays. 

1 Exclusive Nonleptonic B Decays and Factorization 

in : 

\Q . Exclusive nonleptonic decays represent a great challenge to theory. They are 

' complicated by the hadronization of final states and strong-interaction effects 

between them. Today measurements have already reached sufficient precision 
to examine our knowledge of these effects. In order to make real use of data 
. in the determination of fundamental parameters and in testing of the Standard 

Model, we are forced to provide a more accurate estimation of nonperturbative 
quantities, such as the matrix elements of weak operators. 

At the first sight, the nonleptonic B meson decay seems to be simple, as far 
as we essentially consider this decay as a weak decay of heavy b quark. We are 
encouraged to use this argument by the facts that the b quark mass is heavy 
compared to the intrinsic scale of strong interactions and that the b quark decays 
fast enough to produce energetic constituents, which separate without interfer- 
ing with each other. This naive picture was supported by the color-transparency 
argument |EJ and natural application to nonleptonic two-body decays emerged 
under the name the naive factorization (discussed in detail below). However, al- 
though predictions from the naive factorization are in relatively good agreement 
with the data (apart from the color-suppressed decays), the naive factorization 
provides no insight into the dynamical background of exclusive nonleptonic de- 
cays. 

The theoretical discussion of the nonleptonic decay starts with the effective 
weak Hamiltonian, which summarizes our knowledge of weak decays at low scales 
(for a review see 0j): 

H weak = -j= V Qqi V q2q3 [d {jjl)Oi + C 2 (n)<D 2 + ] . (1) 
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The Vs represent the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 
specified for the particular heavy-quark decay Q — > 919293 • Strong-interaction 
effects above some scale /i ~ mb are retained in the Wilson coefficients Ci(n). 
These coefficients are perturbatively calculable and therefore well known. Ac- 
tually, the weak theory without strong corrections and QED effects knows only 
the operator 0i, and in that case C\(Mw) — 1 and Ci{Mw) = 0. The operator 
02, defined in (||), emerges after taking the gluon exchange into account and 
therefore its contribution is suppressed as C 2 (/i) ~ \n(M w ) / \n(fi) . 

The main problem persists in the calculation of matrix elements of operators 

01 in a particular process. In ([!]) we retain only the leading operators 0\ and 2 
and suppress explicitly so called penguin operators, 0^=3, ...,io> Being multiplied 
by, in principle, small Wilson coefficients, the penguin operators usually can be 
neglected (except for the penguin-dominated decays), but could be extremely 
important for detection oi CP violation in B decay 

The four-quark operators 0i and 02 differ only in their color structure: 

Oi = (qu^QiMj^qv) , 02 = (qu^QjKq^qsi) , (2) 

where i and j are color indices, and -T M = 7^(1 — 75). The color-mismatched 
operator 02 can be projected to the color singlet state by using the relation 

SijSki = 1/N c 8«8 jk + 2 {X a /2) tl (X a /2) ]k , as 

02 = ^01+20!. (3) 

This projection, as can be seen from (|^), results in a relative suppression of the 

02 operator contribution of the order 1/N C (N c is the number of colors) and in 
the appearance of the new operator 0\ with the explicit color SU(3) matrices 
A a : 

\ a \ a 

Depending on the process involved, the operators 0\ and 02 can exchange their 
roles, and then it is customary to define the effective parameters a± and ai as 



Ox - QiXq^—qaj) ■ (4) 



ai = d(/x) + i-C 2 (/i), a 2 = C 2 (/i)+^C 1 (/i). (5) 

These parameters distinguish between three classes of decay topologies: 

- class-1 decay amplitude, where a charged meson is directly produced in the 
weak vertex; i.e. in the quark transition b — * udu with 0\ — (dr^^itir^b): 

A(B -> 7T+7T-) ~ ai(Oi) , (6) 

- class-2 decay amplitude, where a neutral meson is directly produced, i.e. in 
the quark transition b — * esc with 2 = (cr^fir^b): 

A(B+ -+ J/ijK+)~a 2 {0 2 ), (7) 
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- class-3 decay amplitude, where both cases are possible, but this amplitude is 
however connected by isospin symmetry with the class-1 and class-2 decays; i.e. 
in the quark transition b — ► csu with 0\ = (ET) J u)(sP''6): 

A{B~ ^D K-)~( ai +xa2){Oi), (8) 

where x denotes the nonperturbative factor being equal to one in the flavor- 
symmetry limit. 

The effective parameters a\ and 02 are defined with respect to the naive 
factorization hypothesis, which assumes that the nonleptonic amplitude can be 
expressed as the product of matrix elements of two hadronic (bilinear) currents, 
for example: 

(n + ir-\(dr^u)(ur»b)\B) -» (tt" | (d/^w) |o) <tt+ | (nr^b) | (9) 

and that there is no nonfactorizable exchange of gluons between the n~ and the 
\ir + B) system. Effectively, that means that the 'nonfactorizable' matrix element 
of the 0\ operator (pj), is vanishing, due to the projection of the colored current 
to the physical colorless state. 



1.1 Nonfactorizable Contributions 

The effective parameters a\ and a-i could be generalized to parametrize also the 
nonfactorizable strong-interaction effects, for example gluon exchanges between 
bilinear currents (i.e. in (^)) which introduce nonvanishing contribution from O 
operators. Schematically, in the large iV c limit, 

01 - Ci(M) + ^(aO + 2<7 2 ( M )£ / (/i) , 

« 2 = C 2 (/i) + ^fM + 2Ci0i)^(/i) , (10) 

where we have explicitly indicated that the nonfactorizable contribution to the 
class-1 and class-2 decays £™l x 2 , do not necessarily need to be the same, and also 
they can be process dependent quantities, which will be discussed later. Theoret- 
ically, nonfactorizable effects are desirable in order to cancel explicit the /i depen- 
dence of Ci(/i) and therefore of the <2j's. All physical quantities are /1 independent, 
and because there is no explicit /i dependence of the matrix elements (Oi) multi- 
plying Ci (fx) , there must be some underlying mechanism to cancel the explicit fi 
dependence of Oj's persisting in the factorization approach. In the calculation of 
the Wilson coefficients beyond the leading order, also the renormalization scheme 
dependence is presented || . Naturally, the parameter a-i is more sensitive on the 
value of the factorization scale and on the renormalization scheme, due to the 
similar magnitude and different sign of the Ci(\i) and \/N c Ci(^jl) terms (calcu- 
lated in the NDR scheme and for A^L = 225 GeV, the Wilson coefficients have 

MS „ 

the following values: C\(mb) = 1.082 and C2(mb) = —0.185 |2j]). This means 
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also that a 2 is more sensitive to any additional nonperturbative long-distance 
contributions. 

The global fit of a\ and a 2 parameters to the B meson experimental data 
performed in |t]], has shown that the a\ coefficient, being essentially proportional 
to Ci(^) ~ 1, is in the expected theoretical range: 

ai ~ 1.05 ±0.10, (11) 

while a 2 has the fitted value of 

a 2 ~ 0.25 ±0.05. (12) 

Compared with the theoretical values calculated with the C\ and C 2 stated 
above, we note that both fitted values show no explicit indication that there is 
a significant nonfactorizable contribution in B decays. This confirms the naive 
factorization picture, although the simple extrapolation of results in D decays 
to the B case would suggest that the 02 coefficient could be negative, meaning 
a nontrivial cancellation of the 1/N C terms and dominance of (negative) C 2 (n) 
in (0). The negative value of a 2 in D decays has found its confirmation in the 
large N c hypothesis of neglecting the higher order 1/N C terms, and in the 
QCD sum rule calculation where the cancellation of the l/N c part with the 
explicitly calculated nonfactorizable terms was verified. 

However, there are additional indications that nonfactorizable contributions 
in B decays cannot be simply neglected and deserve to be investigated. New 
experimental data on B mesons indicate nonuniversality of the a 2 parameter and 
the strong final-state interaction phases in the color-suppressed class-2 decays 
being proportional to a 2 Jl(| ■ 

Therefore, the nonfactorizable contributions must play an important role in 
nonlcptonic decays, particularly in the color-suppressed class-2 decays, such as 
the B — > J/%pK decay discussed in Sect. 4. 

1.2 Models for the Calculation of Nonfactorizable Contributions 

Nowadays, there exist several approaches for the treatment of nonleptonic de- 
cays, which try to investigate the dynamical background and nonfactorizable 
contributions of such processes. The most exploited ones are the QCD-improved 
factorization, fTi"|| , and the PQCD approach []l2"|| . 

The PQCD approach claims the perturbativity of the two-body nonfactoriz- 
able amplitude if the Sudakov suppression is implemented into the calculation. 
The Sudakov form factor suppresses the configuration in which the soft gluon 
exchange could take place, and the amplitude is dominated by exchange of hard 
gluons and therefore perturbatively calculable. 

A somewhat different method is applied in the QCD-improved factorization. 
This method provides the factorization formula that separates soft and hard 
contribution on the basis of large m& expansion. The leading nonfactorizable 
strong interaction can be then studied systematically, while soft (incalculable) 



Exclusive Nonleptonic B decays from QCD Light-Cone Sum Rules 



5 



contributions are suppressed by Aqcd /mb- The method applies to class-1 decays 
and to class-2 decays under the assumption m c <C mi,. 

None of these models take nonfactorizable soft 0(Aqcd /fnb) corrections into 
account. These corrections can be brought under control by using the light-cone 
QCD sum rule method This method is going to be discussed more in detail 
in what follows. 



2 Light-Cone Sum Rules 

All QCD sum rules are based on the general idea of calculating a relevant quark- 
current correlation function and relating it to the hadronic parameters of interest 
via a dispersion relation. Sum rules in hadron physics were already known before 
QCD was established (for a comprehensive introduction to sum rules see i.e. 
[fl4"|), but have reached wide application in a calculation of various hadronic 
quantities in the form of so-called SVZ sum rules fl5|| . The other type of sum 
rules, the light-cone QCD sum rules were established for calculation of exclusive 
amplitudes and form factors (jl6| and references therein). 

2.1 Light-Cone Sum Rules vs SVZ Sum Rules 

In order to illustrate application of the QCD sum rules and the main differences 
between SVZ sum rules and light-cone sum rules, we introduce an example. 

One of the typical calculation using the SVZ sum rules is the estimation of 
the B meson decay constant Jb- The starting point is a correlation function 
defined as 

F(q 2 )=i J d i xe lqx (0\T{7n b mj 5 b(x) 7 m b bi-f 5 u{0)}\0) . (13) 

In the Euclidean region of q momenta, q 2 < 0, we can perform a perturba- 
tive calculation in terms of quarks and gluons by applying the short-distance 
operator-product expansion (OPE) to the correlation function F(q 2 ). The corre- 
lation function is then expressed via a dispersion relation in terms of the spectral 
function p OPE , representing the perturbative part, and the quark and gluon con- 
densates, i.e. (qq), (GG), etc (see for example [fl8||): 

pOPE^ ) = i r ds ImF° PE (s)d S +Ms){omQ) (w) 
n Jo s - T 

where 

-ImF OPE (s)=p OPE (s) (15) 

7T 

and Ai are perturbative coefficients in front of the vacuum condensates of oper- 
ators fli = qq, GG , qX a /2a ■ Gq , etc. 
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On the other hand, in the physical (Minkowskian) region, q 2 > 0, we insert 
the complete sum over hadronic states starting from the ground state B me- 
son, and use a defining relation for fg 
function F(q 2 ) can then be written as 



mi,qi^b\B) = fBm 2 B . The correlation 



F 



hadron 



(9 2 ) 



4 c 
m BJB 



^hadron 



(s) ds 



r 



(16) 



where the hadronic spectral density p hadron contains all higher resonances and 
non-resonant states with the B meson quantum numbers. 

By applying the quark-hadron duality to these higher hadronic (continuum) 
states, which means assuming that we can replace the continuum of hadronic 
states, described by the hadronic spectral function p hadron ( s ^ v ia a dispersion 
relation, by the spectral function calculated perturbatively in the q 2 < re- 
gion p OPE {s), we match both sides, F OPE (q 2 ) = F hadron (q 2 ), and extract the 
needed quantity Jb- The replacement is done for s > s B , where s B is an effec- 
tive parameter of the order of the mass of the first excited B meson resonance 
squared. 

In a practical calculation one performs a finite power expansion in p (s). 
To improve the convergence of the expansion, the Borel transform of both sides, 
F OPE {q 2 ) and F hadron (q 2 ^ [ s considered, defined by the following limiting pro- 
cedure 



B — lira 



1 (- Q ^( — )^ 

q ' [ dq 2> 



— = M fixed . 



(17) 



M 2 is so called Borel parameter. It is determined by the search for stability 
criteria in a sense that, on the one hand, excited and continuum states are sup- 
pressed (asks for smaller M 2 ) and, on the other hand, the reliable perturbative 
calculation is enabled (asks for larger M 2 ). 

The general procedure of QCD sum rules is depicted on Fig. |l|. 



QUARK-HADRON DUALITY 



"quark-gluon side" (f < 



"hadronic side" q">0 



ground 
slate 



ccmliiHium 




Fig. 1. An illustration of the matching procedure in QCD sum rules 
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For calculating quantities which involve hadron interactions, such as for ex- 
ample the B — > 7r form factor, the light-cone sum rules are more suitable ]l7| |. 
The correlation function is now defined as a vacuum-to-pion matrix element: 

= i J d 4 xe-* p ^n(q)\T{u^b(x),m b bi l5 d{0)}\0) . (18) 

The calculation follows by performing a light-cone OPE, an expansion in terms 
of the light-cone wave functions of increasing twist (twist = dimension - spin) . 
Physically, it means that one performs an expansion in the transverse quark dis- 
tances in the infinite momentum frame, rather than a short-distance expansion 
| fl7| . Instead of dealing with the vacuum-to- vacuum quark and gluon condensates 
(numbers) like in the SVZ sum rules, we have now to know the pion distribution 
amplitude (wave function). The leading twist-2 pion distribution amplitude, ^ 
is defined as 

(Tr(q)\u(x)w 5 d(0)\0) = -iq„U I due iu ^^{u). (19) 

Jo 

Distribution amplitudes (DAs) describe distributions of the pion momentum over 
the pion constituents and u denotes the fraction of this momentum, < u < 1 
(for a comprehensive paper on the exclusive decays and the light-cone DAs see 
[fl9|). The DAs represent a nonperturbative, noncalculable input and their form 
has to be determined by nonperturbative methods and/or somehow extracted 
from the experiment. 

In the physical region of (p — q) 2 > nothing changes in comparison to the 
SVZ sum rules. We insert the complete set of hadronic states with B meson 
quantum numbers as before, and extract the B — > tt form factor from the rela- 
tion: {■n(q)\u"f^b\B{p + q)) = 2fg JI (p 2 )q^ l + .... The matching procedure follows 
as described above. 



3 Nonfactorizable Effects in the Light-Cone Sum Rules 

Although the idea to apply QCD sum rules for calculating nonfactorizable con- 
tributions in nonleptonic B decays is not the new one, earlier applications were 
facing some problems which have caused unavoidable theoretical uncertainties 
in their results |l3| ] . In the work |l3| ] , a new approach was introduced and we are 
going first to review its main ideas in the application to the B — > irir decay. 



3.1 Definitions 

The correlator for the B — ► irir decay given in terms of two interpolating currents 
for the pion and the B meson, — w~f v ^d and = mibi^d respectively, 
and relevant operators O x = (dr^ur^b) and O x = (dr^X" '/2u){ur^\ a /2b) 
looks like: 

F„(p,q t k) = J d 4 xe-^-"^ J d i ye^y(0\T{4l\y)O t (0)4 B) (x)}\7T(q)). 

(20) 
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Fig. 2. Sum rule calculation of the B — » tttv decay. The shaded oval region denotes 
nonperturbative input, ir meson distribution amplitude. The other pion and the B 
meson are represented by the currents J^\p — k) and J^ B \p — q) respectively. The 
square represents the four-quark operators Oi. 



The transition is defined again between a vacuum and an external pion state. 
The situation is illustrated in Fig. g. One can note an unphysical momentum 
k coming out from the weak vertex. It was introduced in order to avoid the B 
meson four- momenta before (pb = {p — <?)), and after (P) the decay to be the 
same, Fig. [| In such a way, it was prevented that the continuum of light states 
enters the dispersion relation of the B channel. States, like DD S and D*D S , have 
masses smaller than the ground state B meson mass and spoil the extraction of 
the physical B meson. These 'parasitic' contributions have caused problems in 
the earlier application of the sum rules Jl5| ] . There are several other momenta 
involved into the decay and we take p 2 = k 2 = q 2 = and consider region of 
large spacelike momenta 

\(p-k) 2 \~\(p- q ) 2 \~\P 2 \^A 2 QCD , (21) 

where the correlation function is explicitly calculable. 



3.2 Procedure 

The procedure which one performs is exhibited in Fig. 3. First, Fig. 3a, one 
makes a dispersion relation in a pion channel of momentum (j>— k) 2 and applies 
the quark-hadron duality for this channel, as it was explained in Sect. 2. There- 
after, to be able later to extract physical B meson state, one has to perform an 
analytical continuation of P momentum to its positive value, P 2 = m 2 B . This 
procedure is analogous to the one in the transition from the spacelike to the 
timelike pion form factor, Fig. 3b. Finally, Fig. 3c, a dispersion relation in the 
B channel of momentum (p — q) 2 has to be done, together with the application 
of the quark-hadron duality, now in the B channel ]l3]| . In such a way we ar- 
rive to the double dispersion relation. Apart from somewhat more complicated 
matching procedure, the calculation otherwise follows in a standard way. 



3.3 Results and Implications in the B — > tztz decay 



In 1 1 3 , first, the factorization of the 0\ operator contribution in the B — > 7T7t 
decay was confirmed. The soft nonfactorizable contributions due to the Oi oper- 
ator, which express the exchange of soft gluons between two pions in Fig. @ were 
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(a) Dispersion relation in the pion channel of momentum (p — k) 2 < 
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(b) Analytical continuation of P to P = m B 




JI(-q) 



jt(p) 
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(c) Dispersion relation in the B meson channel of momentum (p — q) 2 < 
Fig. 3. The light-cone sum rule procedure for exclusive decays 
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then calculated. Nonfactorizable soft contributions appear from the absorption 
of a soft gluon emerging from the light-quark loop ud in Fig. ^], by the distri- 
bution amplitude of the outcoming pion 7r + (— q) and there are of the higher, 
twist-3 and twist-4 order in comparison to the factorizable contributions. 

Nonfactorizable soft corrections appeared to be numerically small (~ 1%) 
and suppressed by l/m^. Therefore, their impact on the complete decay ampli- 
tude was shown to be of the same order as that of hard nonfactorizable contri- 
butions calculated in the QCD-improved factorization approach fll| . Also, the 
calculation has shown no imaginary phase from the soft contributions, whereas 
aforementioned hard nonfactorizable contributions get small complex phase be- 
cause of the final state rescattering due to the hard gluon exchange. 



4 Nonfactorizable Effects for B — > J/iftK 

The B — > J/tpK decay was considered in |2(J. As it was emphasized at the 
beginning, this decay belongs to the color-suppressed class-2 decays in which one 
expects large nonfactorizable contributions. The confirmation of this assumption 
seems to be also found experimentally. Namely, there is a discrepancy between 
the experiment and the naive factorization prediction by at least a factor of 3 in 
the branching ratio. The Hamiltonian which describes the decay is given as 



Gf 

H w = —j=V c bV* t 



(C 2 ( A1 ) + ^C 1 ( Ai ))0 2 + 2C 1 ( A1 )a 2 



(22) 



with the operators 2 = (cr^isr^b) and 6 2 = (aT,, 4fc)(s/^4f 6). In the 
factorization approach, the matrix element of Oi vanishes, and the factorized 
matrix element of the operator O2 is given by 

(J/i>(p)K(q)\0 2 \B(p + q)) = (J/^(p)\cr^\0){K(q)\sr»b\B(p + q)) 

= 2e-q mj/^fj^F+^m^j^) . (23) 

Fg K {m 2 j^) is the B — ► K transition form factor calculated using the light-cone 
sum rules, in a way enlightened in Sect. 2.1 on the example of B — > n form factor 
calculation, and /ju is the J/ip decay constant. By evaluating numerically the 
B — > J/tpK branching ratio with the NLO Wilson coefficients used in Sec. 1.1 
and with the numerical input taken from p0[ , we arrive to 

B(B -» Jji\)K) Sact = 3.3 • 1(T 4 . (24) 

This has to be compared with the recent measurements |2l| 

B{B+ -> J/ipK + ) = (10.1 ± 0.3 ± 0.5) • 10~ 4 , 

B(B° -> J/^jK ) = (8.3 ±0.4 ±0.5) • 10~ 4 . (25) 



It is clear that there a discrepancy between the naive factorization prediction 
and the experiment. 



Exclusive Nonleptonic B decays from QCD Light-Cone Sum Rules 11 

To be able to discuss the impact of the nonfactorizable term 2 , we parametrize 
the ( J/ipK\Hw\B) amplitude in terms of the a 2 parameter as 

(J/iPK[H w \B) =V2G F V cb V c * s e • qmj^fj^F+^m 2 ^) a 2 , (26) 

where 

a 2 = C 2 ( P ) + ^ + 2d ( M ) : • (27) 

The part proportional to Fg K represents the contribution from the <D 2 operator 
(J/iPK\6 2 (p)\B) = 2e ■ 9 m w / J/ /+ Jf ^) (28) 



and Fg K = corresponds to the naive factorization result, Eq. (p3[). 

By using the parametrization (|2^) we can extract the a 2 coefficient from 
experiments (|25|). The measurements yield 

\af p \ =0.29 ±0.03. (29) 

On the other hand, the naive factorization with the NLO Wilson coefficients 
[|J produces 

fl2°«LO = °' 176 ■ (30) 

The value (|30|) is significantly below the value extracted from the experiment, 
although one should not forget a strong /z dependence of a 2 act . 

In Fig. [| we show the partial width for B — > J/ipK as a function of the 
nonfact 




Fig. 4. The partial width for B — > J/ipK as a function of the nonfactorizable amplitude 
Fsk- The dashed-dotted lines denote the experimental region. 
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factorizable prediction. There exist two ways to satisfy the experimental de- 
mands on Fbk- Following the large 1/N C rule Q, one can argue that there is a 
cancellation between 1 /N c piece of the factorizable part and the nonfactorizable 



contribution (27). This would ask for the relatively small and negative value of 
Fbk- The other possibility is to have even smaller, but positive values for Fbk, 
which then compensate the overall smallness of the factorizable part and bring 
the theoretical estimation for a?, in accordance with experiment. 

One can note significant [i dependence of the theoretical expectation for the 
partial width in Fig. [|, which brings an uncertainty in the prediction for Fbk{h) 
in the order of 30%. This uncertainty is even more pronounced for the positive 
solutions of Fbk (m) ■ The values for Fg K extracted from experiments 

F+ K (m b ) = 0.028 or F+ K (m b ) = -0.120 , (31) 
F+ K (m b /2) = 0.046 or F^ K (m b /2) = -0.095 . (32) 

clearly illustrate the /i sensitivity of the nonfactorizable part. 

In what follows we calculate the nonfactorizable contribution F^ K which 
appears due to the exchange of soft gluons using the QCD light-cone sum rule 
method. 



4.1 Light-Cone Sum Rule Calculation 

The light-cone sum rule calculation starts by considering the correlator 

F lJ (p,q,k)=2 2 J d 4 xe~^ + ^ J d%e^-Vy(K(q)\T{jy^Xy)O(0)4 B \x)}\0) 

(33) 

with the interpolating currents Jv 1 ^' > = cj„c and jj- B ^ = m b bi^u. The kinemat- 
ics is the same as defined above in (^lj), with the exception that now p 2 = m 2 ^. 
More explicitly the configuration is shown in Fig. 

The estimation of nonfactorizable contributions was performed for the ex- 
change of soft gluons (shown by the dashed line in Fig. ||) and follows essentially 

K(q) 




Fig. 5. Sum rule calculation of the B — > J/ipK decay. The dashed line denote an 
exchange of a soft gluon for Oi = O2 and the cross stands for other possible attachment 
of a soft gluon. 
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steps of derivation explained in Sect. 3.2 for B — ► irir decay. Nonfactorizable 
contribution of the 0\ operator appears first at 0{a 2 s ). Nonvanishing result at 
the one-gluon level includes contribution of the O2 operator and the leading con- 
tributions are given in terms of twist-3 and twist-4 kaon distribution amplitudes 
which contribute in the same order. Technical peculiarities of the calculation can 
be found in 

4.2 Results and Implications 

The results can be summarized as follows. Soft nonfactorizable twist-3 and twist- 
4 contributions, expressed in terms of Fg K are Fg K tw3 (fJ-b) = 0.003 — 0.0055 
and Fg K twi {^b) = 0.006 — 0.0012 and the final value is 

f bkM = 0.009 - 0.0017 . (34) 

where pn, — \J Mg — rri^ ~ TOb/2. The wide range prediction for Fg K appears 
due to the variation of sum rule parameters. 

First, we note that the nonfactorizable contribution ( |34|) is much smaller 
than the B — > K transition form factor Fg K — 0.55 ± 0.05, which enters the 
factorization prediction (26). It is also significantly smaller than the value ([32J) 
extracted from experiments. Nevertheless, its influence on the final prediction 
for 02 is significant, because of the large coefficient 2C\ multiplying it. Further, 
one has to emphasize that Fi K is a positive quantity. Therefore, we do not find 
a theoretical support for the large iV c limit assumption discussed in Sect. 4.1, 
that the factorizable part proportional to Ci(/i)/3 should at least be partially 
cancelled by the nonfactorizable part. Our result also contradicts the result of 
the earlier application of QCD sum rules to B — > J/t/jK S), where negative 
and somewhat larger value for Fg K was found. However, earlier applications of 
QCD sum rules to exclusive B decays exhibit some deficiencies discussed in ]l3]| . 

Using the same values for the NLO Wilson coefficients as in Sect. 2, one gets 
from (|34|) the following value for the effective coefficient 02 : 

a 2 ~ 0.14- 0.17 | M=W . (35) 

Although the soft correction contributes in the order of ~ 30% — 70%, the net 
result ([[5]) is still by approximately factor of two smaller than the experimentally 
determined value (|29|). 

We would like to discuss our results for soft nonfactorizable contributions 
in comparison with the hard nonfactorizable effects calculated in the QCD- 
improved factorization approach. The best thing would be to calculate both soft 
and hard contributions inside the same model. In principle, the light-cone sum 
rule approach presented here enables such a calculation, although the estimation 
of hard nonfactorizable contributions is technically very demanding, involving a 
calculation of two-loop diagrams. Therefore, we proceed with the QCD-improved 
factorization estimations for the hard nonfactorizable contributions. 
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After including the hard nonfactorizable corrections, the 02 parameter (|35| 
is as follows 



a-2 



BK 



(36) 



The estimations done in the QCD-improved factorization ]24| show hard- 
gluon exchange corrections to the naive factorization result in the order of ~ 25%, 
predicted by the LO calculation with the twist-2 kaon distribution amplitude. 
Unlikely large corrections are obtained by the inclusion of the twist-3 kaon distri- 
bution amplitude. Anyhow, due to the obvious dominance of soft contributions 
to the twist-3 part of the hard corrections in the BBNS approach |Tj|| , it is 
very likely that some double counting of soft effects could appear if we naively 
compare the results. Therefore, taking only the twist-2 hard nonfactorizable cor- 
rections from Q into account, recalculated at the fib scale, our prediction ( ^5| ) 
changes to 

02 = 0.16 - 0.19 | p=w (37) 

The prediction still remains too small to explain the data. 

Nevertheless, there are several things which have to be stressed here in con- 
nection with the result. Soft nonfactorizable contributions are at least equally 
important as nonfactorizable contributions from the hard-gluon exchange, and 
can be even dominant. Soft nonfactorizable contributions are positive, and the 
same seems to be valid for hard corrections. While hard corrections have an imag- 
inary part, in the soft contributions the annihilation and the penguin topologies 
as potential sources for the appearance of an imaginary part were not discussed. 
A comparison between the result (EjTl) and the experimental value I02I ~ 0.3 
for B — > J /ipK decay, with the recently deduced 02 parameter from the color- 
suppressed B — ► _DW°7r° decays, |oa| ~ 0.4 — 0.5 [|l0), provides clear evidence 
for the nonuniversality of the ai parameter in color-suppressed decays. 



5 Conclusions 

We have reviewed recent progress in the understanding of the underlaying dy- 
namics of exclusive nonleptonic decays, with the emphasis on the nonfactorizable 
corrections to the naive factorization approach. In the calculation of nonfactor- 
izable contributions, we have focused to QCD light-cone sum rule approach and 
have shown results for B — > -kit [[l3| and B — > J/ipK [^0) decays. 

The QCD-improved factorization method is reviewed in this volume by M. 
Neubert, j§. 
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Abstract. We are going to review recent advances in the theory of exclusive nonlep- 
tonic B decays. The emphasis is going to be on the factorization hypothesis and the role 
of nonfactorizable contributions for nonleptonic B decays. In particular, we will discuss 
in more detail calculations of nonfactorizable contributions in the QCD light-cone sum 
rule approach and their implications to the B — > tttt and B —¥ J/ipK decays. 

1 Exclusive Nonleptonic B Decays and Factorization 

Exclusive nonleptonic decays represent a great challenge to theory. They are 
complicated by the hadronization of final states and strong-interaction effects 
between them. Today measurements have already reached sufficient precision 
to examine our knowledge of these effects. In order to make real use of data 
in the determination of fundamental parameters and in testing of the Standard 
Model, we are forced to provide a more accurate estimation of nonperturbative 
quantities, such as the matrix elements of weak operators. 

At the first sight, the nonleptonic B meson decay seems to be simple, as far 
as we essentially consider this decay as a weak decay of heavy b quark. We are 
encouraged to use this argument by the facts that the b quark mass is heavy 
compared to the intrinsic scale of strong interactions and that the b quark decays 
fast enough to produce energetic constituents, which separate without interfer- 
ing with each other. This naive picture was supported by the color-transparency 
argument [1] and natural application to nonleptonic two-body decays emerged 
under the name the naive factorization (discussed in detail below). However, al- 
though predictions from the naive factorization are in relatively good agreement 
with the data (apart from the color-suppressed decays), the naive factorization 
provides no insight into the dynamical background of exclusive nonleptonic de- 
cays. 

The theoretical discussion of the nonleptonic decay starts with the effective 
weak Hamiltonian, which summarizes our knowledge of weak decays at low scales 
(for a review see [2]): 

n weak = % v Qqi v q2qs [Ci(a*)Oi + c 2 (n)o 2 + ] . (i) 
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The Vs represent the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 
specified for the particular heavy-quark decay Q — > qiq 2 q 3 . Strong-interaction 
effects above some scale [i ~ to& are retained in the Wilson coefficients Ci(fi). 
These coefficients are perturbatively calculable and therefore well known. Ac- 
tually, the weak theory without strong corrections and QED effects knows only 
the operator 0i, and in that case C\(Mw) = 1 and C 2 (Mw) = 0. The operator 
02, defined in (2), emerges after taking the gluon exchange into account and 
therefore its contribution is suppressed as C 2 (fi) ~ ln(Mw)/ ln(/i). 

The main problem persists in the calculation of matrix elements of operators 

01 in a particular process. In (1) we retain only the leading operators 0\ and 02 
and suppress explicitly so called penguin operators, 0j=3,...,io. Being multiplied 
by, in principle, small Wilson coefficients, the penguin operators usually can be 
neglected (except for the penguin-dominated decays), but could be extremely 
important for detection of CP violation in B decay [3-5] . 

The four-quark operators 0\ and 02 differ only in their color structure: 

Oi = (q li r„Q i )(q 2j r'* q3j ) , 02 = (q li r ll Q j )fa j r»q 3i ) , (2) 

where i and j are color indices, and = 7 M (1 — 75). The color-mismatched 
operator 02 can be projected to the color singlet state by using the relation 
SijSu = l/N c S u S jk + 2 (\ a /2) it (\ a /2) jk , as 

02 = ^01+20!. (3) 

This projection, as can be seen from (3), results in a relative suppression of the 

02 operator contribution of the order 1/N C (N c is the number of colors) and a 
new operator 0\ with the explicit color SU(3) matrices A°: 

Oi = (q u r,Y^^ r "Y q3j) - (4) 

Depending on the process involved, the operators 0\ and 02 can exchange their 
roles, and then it is customary to define the effective parameters a\ and ci2 as 

ai = CiM + jJrC 2 (p), a 2 = C 2 ((i) + j^dhi) . (5) 

These parameters distinguish between three classes of decay topologies: 

- class-1 decay amplitude, where a charged meson is directly produced in the 
weak vertex; i.e in the quark transition b — > udu with 0\ = (dr^^r^b): 

A(B ^tt+tt") ~ai<Oi>, (6) 

- class-2 decay amplitude, where a neutral meson is directly produced, i.e in the 
quark transition b — > esc with 02 = (cr^^sT^b): 



A(B+ ^J/iPK+)~a 2 (0 2 ), 



(7) 
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- class-3 decay amplitude, where the both cases are possible, but this amplitude 
is however connected by isospin symmetry with the class-1 and class-2 decays; 
i.e in the quark transition b — > csu with 0\ = (cr^^sT^b): 

A{B- D°K-) ~ (ai + ia 2 )<Oi> , (8) 

where x denotes the nonperturbative factor being equal to one in the flavor- 
symmetry limit. 

The effective parameters a\ and a 2 are defined with respect to the naive 
factorization hypothesis, which assumes that the nonleptonic amplitude can be 
expressed as the product of matrix elements of two hadronic (bilinear) currents, 
for example: 

<7r + 7r-|(dJ>)(«r"&)|fl) -> <7T" | (d/^u) 1 0> <7r+ 1 («^6) I (9) 

and that there is no nonfactorizable exchange of gluons between the 7r~ and the 
\ir + B) system. Effectively, that means that the 'nonfactorizable' matrix element 
of the 0\ operator (4), is vanishing, due to the projection of the colored current 
to the physical colorless state. 

1.1 Nonfactorizable Contributions 

The effective parameters a\ and a 2 could be generalized to parametrize also the 
nonfactorizable strong-interaction effects, for example gluon exchanges between 
bilinear currents (i.e. in (9)) which introduce nonvanishing contribution from O 
operators. Schematically, in the large N c limit, 

a, = KM + ^-C 3 (fi) + 2C 2 (p)&(p) , 

a 2 = C 2 (p) + ^-Ci(fi) + 2C 1 (/i)£ / (/i) , (10) 

where we have explicitly indicated that the nonfactorizable contribution to the 
class-1 and class-2 decays 2 , do not necessarily need to be the same, and also 
they can be process dependent quantities, which will be discussed later. Theoret- 
ically, nonfactorizable effects are desirable in order to cancel explicit the /x depen- 
dence of Ci(fi) and therefore of the a,'s. All physical quantities are [i independent, 
and because there is no explicit fi dependence of the matrix elements (0,) multi- 
plying Ci (p) , there must be some underlying mechanism to cancel the explicit [i 
dependence of a,'s persisting in the factorization approach. In the calculation of 
the Wilson coefficients beyond the leading order, also the renormalization scheme 
dependence is presented [6]. Naturally, the parameter a 2 is more sensitive on the 
value of the factorization scale and on the renormalization scheme, due to the 
similar magnitude and different sign of the C 2 {n) and l/7V c Ci(/i) terms (calcu- 
lated in the NDR scheme and for A^— = 225 GeV, the Wilson coefficients have 

MS ' 

the following values: Ci(m&) = 1.082 and C 2 (mi,) = —0.185 [2]). This means 
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also that a-i is more sensitive to any additional nonperturbative long-distance 
contributions. 

The global fit of a\ and a-i parameters to the B meson experimental data 
performed in [7], has shown that the a\ coefficient, being essentially proportional 
to Ci(/i) ~ 1, is in the expected theoretical range: 

en ~ 1.05 ±0.10, (11) 

while ai has the fitted value of 

a 2 ~ 0.25 ±0.05. (12) 

Compared with the theoretical values calculated with the C\ and C2 stated 
above, we note that both fitted values show no explicit indication that there is 
a significant nonfactorizable contribution in B decays. This confirms the naive 
factorization picture, although the simple extrapolation of results in D decays 
to the B case would suggest that the ai coefficient could be negative, meaning 
a nontrivial cancellation of the 1 /N c terms and dominance of (negative) C2 {p) 
in (10). The negative value of in D decays has found its confirmation in the 
large N c hypothesis of neglecting the higher order 1/N C terms, [8], and in the 
QCD sum rule calculation [9], where the cancellation of the 1/N C part with the 
explicitly calculated nonfactorizable terms was verified. 

However, there are additional indications that nonfactorizable contributions 
in B decays cannot be simply neglected and deserve to be investigated. New 
experimental data on B mesons indicate nonuniversality of the a-i parameter and 
the strong final-state interaction phases in the color-suppressed class-2 decays 
being proportional to ai [10]. 

Therefore, the nonfactorizable contributions must play an important role in 
nonleptonic decays, particularly in the color-suppressed class-2 decays, such as 
the B — » J/tjjK decay discussed in Sect. 4. 

1.2 Models for the Calculation of Nonfactorizable Contributions 

Nowadays, there exist several approaches for the treatment of nonleptonic de- 
cays, which try to investigate the dynamical background and nonfactorizable 
contributions of such processes. The most exploited ones are the QCD-improved 
factorization, [11], and the PQCD approach [12]. 

The PQCD approach claims the perturbativity of the two-body nonfactoriz- 
able amplitude if the Sudakov suppression is implemented into the calculation. 
The Sudakov form factor suppresses the configuration in which the soft gluon 
exchange could take place, and the amplitude is dominated by exchange of hard 
gluons and therefore perturbatively calculable. 

A somewhat different method is applied in the QCD-improved factorization. 
This method provides the factorization formula that separates soft and hard 
contribution on the basis of large mj expansion. The leading nonfactorizable 
strong interaction can be then studied systematically, while soft (incalculable) 
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contributions are suppressed by AQCD/rrib- The method applies to class- 1 decays 
and to class-2 decays under the assumption m c <C mj. 

None of these models take nonfactorizable soft 0(Aqcd /mb) corrections into 
account. This corrections can be brought under the control by using the light- 
cone QCD sum rule method [13]. This method is going to be discussed in more 
details in what follows. 

2 Light-Cone Sum Rules 

All QCD sum rules are based on the general idea of calculating a relevant quark- 
current correlation function and relating it to the hadronic parameters of interest 
via a dispersion relation. Sum rules in hadron physics were already known before 
QCD was established (for a comprehensive introduction to sum rules see i.e. 
[14]), but have reached wide application in a calculation of various hadronic 
quantities in the form of so-called SVZ sum rules [15]. The other type of the 
sum rules, the light-cone QCD sum rules were established for calculation of 
exclusive amplitudes and form factors ([16] and references therein). 

2.1 Light-Cone Sum Rules vs SVZ Sum Rules 

In order to illustrate application of the QCD sum rules and the main differences 
between SVZ sum rules and light-cone sum rules, we introduce an example. 

One of the typical calculation using the SVZ sum rules is the estimation of 
the B meson decay constant The starting point is a correlation function 
defined as 

F(q 2 )=i j d 4 xe iqx (0\T{m b ui-/ 5 b(x),rn b bi-/ 5 u(0)}\0). (13) 

In the Euclidean region of q momenta, q 2 < 0, we can perform a perturba- 
tive calculation in terms of quarks and gluons by applying the short-distance 
operator-product expansion (OPE) to the correlation function F(q 2 ). The corre- 
lation function is then expressed via a dispersion relation in terms of the spectral 
function p OPE , representing the perturbative part, and the quark and gluon con- 
densates, i.e (qq), (GG), etc (see for example [18]): 

F oPE {q 2 ) = ± r ds ImF ° PE ^ ds +M s)(O\fi i \0), (14) 
TT Jo S ~Q 

where 

-ImF OPE (s) = p OPE (s) (15) 

7T 

and Ai are perturbative coefficients in front of the vacuum condensates of oper- 
ators fli = qq, GG , q\ a /2a ■ Gq , etc. 
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On the other hand, in the physical (Minkowskian) region, q 2 > 0, we insert 
the complete sum over hadronic states starting from the ground state B me- 
son, and use a defining relation for fg : (mbqi^5b\B) = fBm 2 B . The correlation 
function F(q 2 ) can be then written as 

4 f2 roc hadron/ „\ J„ 
phadron,^ = J^B^ + / P IfJff (lg) 

m 2 b - q 2 J s h s-q 2 

where the hadronic spectral density p hadron contains all higher resonances and 
non-resonant states with the B meson quantum numbers. 

By applying the quark-hadron duality to these higher hadronic (continuum) 
states, which means assuming that we can replace the continuum of hadronic 
states, described by the hadronic spectral function p hadron ^ v i a a dispersion 
relation, by the spectral function calculated perturbatively in the q 2 < re- 
gion p OPE (s), we match both sides, F OPE (q 2 ) = F hadron (q 2 ), and extract the 
needed quantity The replacement is done for s > Sq , where s B is an effec- 
tive parameter of the order of the mass of the first excited B meson resonance 
squared. 

In a practical calculation one performs a finite power expansion in p OPE (s). 
To improve the convergence of the expansion, the Borel transform of both sides, 
F OPE (q 2 ) and F hadron (q 2 ^ [ s considered, defined by the following limiting pro- 
cedure 

B = lim—^— -(- g 2 )«(A)« \q% n ^ c,^ = M 2 fixed. (17) 
(n — 1)! aq z n 

M 2 is so called Borel parameter. It is determined by the search for stability 
criteria in a sense that, on the one hand, excited and continuum states are sup- 
pressed (asks for smaller M 2 ) and, on the other hand, the reliable perturbative 
calculation is enabled (asks for larger M 2 ). 

The general procedure of QCD sum rules is depicted on Fig. 1. 



QUARK-HADRON DUALITY 




B B' $ 



Fig. 1. An illustration of the matching procedure in QCD sum rules 
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For calculating quantities which involve hadron interactions, such as for ex- 
ample the B — > 7T form factor, the light-cone sum rules are more suitable [17]. 
The correlation function is now defined as a vacuum-to-pion matrix element: 

F M = i j d 4 xe- ip ^7r(q)\T{u^b(x),m b bi 75 d(0)}\0) . (18) 

The calculation follows by performing a light-cone OPE, an expansion in terms 
of the light-cone wave functions of increasing twist (twist = dimension - spin) . 
Physically, that means that one performs an expansion in the transverse quark 
distances in the infinite momentum frame, rather than a short-distance expan- 
sion [17]. Instead dealing with the vacuum-to- vacuum quark and gluon con- 
densates (numbers) as in the SVZ sum rules, we have now to know the pion 
distribution amplitude (wave function). The leading twist-2 pion distribution 
amplitude, <j>„ is defined as 

(7r(q)\u(x) W d(0)\0) = -iq^U [ due iu " x ^{u). (19) 

Jo 

Distribution amplitudes (DAs) describe distributions of the pion momentum over 
the pion constituents and u denotes the fraction of this momentum, < u < 1 
(for a comprehensive paper on the exclusive decays and the light-cone DAs see 
[19]). The DAs represent a nonperturbative, noncalculable input and their form 
has to be determined by nonperturbative methods and/or somehow extracted 
from the experiment. 

In the physical region of (p — q) 2 > nothing changes in comparison to the 
SVZ sum rules. We insert the complete set of hadronic states with B meson 
quantum numbers as before, and extract the B — > 7r form factor from the rela- 
tion: (ir(q)\Ti'y l j / b\B(p + q)) = 2fg 7T (p 2 )q IJ + .... The matching procedure follows 
as described above. 

3 Nonfactorizable Effects in the Light-Cone Sum Rules 

Although the idea to apply QCD sum rules for calculating nonfactorizable con- 
tributions in nonleptonic B decays is not the new one, earlier applications were 
facing some problems which have caused unavoidable theoretical uncertainties 
in their results [13]. In the work [13], a new approach was introduced and we are 
going first to review its main ideas in the application to the B — > tttt decay. 

3.1 Definitions 

The correlator for the B — > tttt decay given in terms of two interpolating currents 
for the pion and B meson, = uj^j^d and = mObi^^d respectively, and 
relevant operators d = (dr^ur^b) and O x = (dr^\ a /2u)(ur»\ a /2b) looks 
like: 

F v {p,q,k) = j d 4 xe-^-^ J d 4 i/e^-*)»<0|T{jW(i/)O i (0)J B (B) (i)}|7r(g)>. 

(20) 
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Fig. 2. Sum rule calculation of the B — > iziz decay. The shaded oval region denotes 
nonperturbative input, 7r meson distribution amplitude. The other pion and the B 
meson are represented by the currents J^'(p — k) and J*- B '(p — q) respectively. The 
square represents the four-quark operators Oi. 



The transition is again defined between a vacuum and an external pion state. 
The situation is illustrated in Fig. 2. One can note an unphysical momentum 
k coming out from the weak vertex. It was introduced in order to avoid the 
B meson four-momenta before (ps = (p — (?)), and after (P) the decay to be 
the same, Fig. 2. In such a way, it was prevented that the continuum of light 
states enters the dispersion relation of the B channel. Contributions, like DD S 
and D*D S , have masses smaller than the ground state B meson mass and spoil 
the extraction of the physical B state. Contributions of this 'parasitic' states 
have caused problems in the earlier application of the sum rules [13]. There are 
several other momenta involved into the decay and we take p 2 = k 2 = q 2 = 
and consider region of large spacelike momenta 

\(p-k) 2 \~\(p-q) 2 \~\P 2 \»A 2 QCD , (21) 

where the correlation function is explicitly calculable. 

3.2 Procedure 

The procedure which one performs is exhibited in Fig. 3.1. First, Fig. 3a, one 
makes a dispersion relation in a pion channel of momentum (p — k) 2 and applies 
the quark-hadron duality for this channel, as it was explained in Sect. 2. There- 
after, to be able later to extract physical B meson state, one has to perform an 
analytical continuation of P momentum to its positive value, P 2 = rri 2 B . This 
procedure is analogous to the one in the transition from the spacelike to the 
timelike pion form factor, Fig. 3b. Finally, Fig. 3c, a dispersion relation in the B 
channel of momentum (p — q) 2 has to be done, together with the application 
of the quark-hadron duality, now in the B channel [13]. In such a way we ar- 
rive to the double dispersion relation. Apart from somewhat more complicated 
matching procedure, the calculation otherwise follows in a standard way. 

3.3 Results and Implications in the B — >• wir decay 

In [13], first the factorization of the 0\ operator contribution in the B — > tttt 
decay was confirmed. The soft nonfactorizable contributions due to the 0\ oper- 
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P = (p-q-kf < o 

(a) Dispersion relation in the pion channel of momentum (p — k) 2 < 




P = (p-k-qj 2 > 



(c) Dispersion relation in the B meson channel of momentum (p — q) 2 < 
Fig. 3. The light-cone sum rule procedure for exclusive decays 
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ator, which express the exchange of soft gluons between two pions in Fig. 2 were 
then calculated. Nonfactorizable soft contributions appear from the absorption 
of a soft gluon emerging from the light-quark loop ud in Fig. 2, by the distri- 
bution amplitude of the outcoming pion 7r + (— q) and there are of the higher, 
twist-3 and twist-4 order in comparison to the factorizable contributions. 

Nonfactorizable soft corrections appeared to be numerically small (~ 1%) 
and suppressed by l/m&. Therefore, their impact to the complete decay ampli- 
tude was shown to be of the same order as that of hard nonfactorizable contri- 
butions calculated in the QCD-improved factorization approach [11]. Also, the 
calculation has shown no imaginary phase from the soft contributions, whereas 
aforementioned hard nonfactorizable contributions get small complex phase be- 
cause of the final state rescattering due to the hard gluon exchange. 



4 Nonfactorizable Effects for B — > J/tjiK 

The B — > J/ipK decay was considered in [20]. As it was emphasized at the 
beginning, this decay belongs to the color-suppressed class-2 decays in which one 
expects large nonfactorizable contributions. The confirmation of this assumption 
seems to be also found experimentally. Namely, there is a discrepancy between 
the experiment and the naive factorization prediction by at least a factor of 3 in 
the branching ratio. The Hamiltonian which describes the decay is given as 



Gf 



(22) 



with the operators 2 = {cr„c){sr^b) and 6 2 = (cr^c)(sr»^b). In the 
factorization approach, the matrix element of 2 vanishes, and the factorized 
matrix element of the operator 2 is given by 

(J/^{p)K(q)\0 2 \B{p + q)) = (J/*[,(p)\cr^c\0)(K(q)\sr»b\B(p + q)) 

= 2e ■qmj /l pfj /l pF+ K (m 2 j /lp ) . (23) 

^sif ( m j/«/>) 1S the- B — > K transition form factor calculated using the light-cone 
sum rules, in a way enlightened in Sect. 2.1 on the example of B — > 7r form 
factor calculation, and fj/^ is J/tjj decay constant. By evaluating numerically 
the B — > J /ipK branching ratio with the NLO Wilson coefficients used in Sec. 1.1 
and with the numerical input taken from [20], we arrive to 

B(B Jl4>K) fact = 3.3 • 10" 4 . (24) 

This has to be compared with the recent measurements [21] 

B(B + JH>K+) = (10.1 ± 0.3 ± 0.5) • 10" 4 , 

B(B° ->■ J/tpK°) = (8.3 ±0.4 ±0.5) • 10" 4 . (25) 



It is clear that there a discrepancy between the naive factorization prediction 
and the experiment. 
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To be able to discuss the impact of the nonfactorizable term O2 , we parametrize 
the (J/ipK \Hw\B) amplitude in terms of the a-i parameter as 

(J/^K\H W \B) =V2G F V cb V c * s e ■ qmj^fj^F+ K (m 2 J/jp ) a 2 , (26) 

where 

« 2 = C 2 M + + 2Cl fo) . (27) 

The part proportional to Fg K represents the contribution from the 2 operator 

{Jl^K\6 2 {n)\B) = 2e ■qmj / ^fj^F+ K {i?) (28) 

and Fg K = corresponds to the naive factorization result, Eq. (23). 

By using the parametrization (26) we can extract the ci2 coefficient from 
experiments (25). The measurements yield 



exp i 



0.29 ± 0.03. (29) 



On the other hand, the naive factorization with the NLO Wilson coefficients 
[6] produces 

°%NLO = °- 176 Usims • ( 30 ) 

The value (30) is significantly below the value extracted from the experiment, 
although one should not forget a strong \i dependence of a$^ ct . 

In Fig. 4 we show the partial width for B — > J/ipK as a function of the 
nonfactorizable amplitude Fb^k- The zero value of Fbk corresponds to the 




Fig. 4. The partial width for B — > J/ipK as a function of the nonfactorizable amplitude 
Fbk- The dashed-dotted lines denote the experimental region. 



12 Blazenka Melic 



factorizable prediction. There exist two ways to satisfy the experimental de- 
mands on Fbk ■ Following the large 1 /N c rule [8] , one can argue that there is a 
cancellation between 1 /N c piece of the factorizable part and the nonfactorizable 
contribution, (27). This would ask for the relatively small and negative value of 
Fbk • The other possibility is to have even smaller, but positive values for Fbk, 
which then compensate the overall smallness of the factorizable part and bring 
the theoretical estimation for cii in agreement with experiment. 

One can note significant /U dependence of the theoretical expectation for the 
partial width in Fig. 4, which brings an uncertainty in the prediction for Fbk{p) 
of the order of 30%. This uncertainty is even more pronounced for the positive 
solutions of Fbk (/■*) • The values for Fg K extracted from experiments 

F+ K (m b ) = 0.028 or F+ K (m b ) = -0.120 , (31) 
F+ K (m b /2) = 0.046 or F+ K (m b ) = -0.095. (32) 

clearly illustrate the /x sensitivity of the nonfactorizable part. 

In the following we will calculate the nonfactorizable contribution Fg K which 
appears due to the exchange of soft gluons using the QCD light-cone sum rule 
method. 



4.1 Light-Cone Sum Rule Calculation 

The light-cone sum rule calculation starts by considering the correlator 

F v (p, q, k)=i 2 J d i xe- i{ P + "> J ^y^-^ (K(q)\T{ji J l^\y)O(0)4 B \x)}\0) 

(33) 

with the interpolating currents jl J ^ = c^ v c and Jg = m b bij^u. The kinemat- 
ics is the same as defined above in (21), with the exception that now p 2 = rn?j^. 
More explicitly the configuration is shown in Fig. 5. 

The estimation of nonfactorizable contributions was performed for the ex- 
change of soft gluons (shown by the dashed line in Fig. 5) and follows essentially 




Fig. 5. Sum rule calculation of the B — > J/ipK decay. The dashed line denote an 
exchange of a soft gluon for Oi = O2 and the cross stands for other possible attachment 
of a soft gluon. 
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steps of derivation explained in Sect. 3.2 for B — > irir decay. Nonfactorizable 
contribution of the 0\ operator appears first at 0(a 2 s ). Nonvanishing result at 
the one-gluon level includes contribution of the O2 operator and the leading con- 
tributions are given in terms of twist-3 and twist-4 kaon distribution amplitudes 
which contribute at the same order. Technical peculiarities of the calculation can 
be found in [20]. 

4.2 Results and Implications 

The results can be summarized as follows. Soft nonfactorizable twist-3 and twist- 
4 contributions, expressed in terms of Fg K are Fg K tw3 (lJ>b) = 0.003 — 0.0055 
and fi£ K tw4 (ii b ) = 0.006 - 0.0012 and the final value is 

F+ K (fi b ) = 0.009 - 0.0017 . (34) 

where /ii, = yj Af| — ~ toj/2. The wide range prediction for Fg K appears 
due to the variation of sum rule parameters. 

First, we note that the nonfactorizable contribution (34) is much smaller 
than the B — > K transition form factor Fg K = 0.55 ± 0.05, which enters the 
factorization prediction (26). It is also significantly smaller than the value (32) 
extracted from experiments. Nevertheless, its influence on the final prediction 
for ai is significant, because of the large coefficient 2C\ multiplying it. Further, 
one has to emphasize that Fg K is a positive quantity. Therefore, we do not 
find a theoretical support for the large N c limit assumption discussed in Sec.4.1, 
that the factorizable part proportional to Ci(/i)/3 should at least be partially 
canceled by the nonfactorizable part. Our result also contradicts the result of 
the earlier application of QCD sum rules to the B — > J/ipK [23], where negative 
and somewhat larger value for Fg K was found. However, earlier applications of 
the QCD sum rules to exclusive B decays exhibit some deficiencies discussed in 
[13]. 

Using the same values for the NLO Wilson coefficients as in Sec. 2, one gets 
from (34) the following value for the effective coefficient a^: 

a 2 -0.14 -0.17 (35) 

Although the soft correction contribute at the order of ~ 30% — 70%, the net 
result (35) is still by approximately factor of two smaller than the experimentally 
determined value (29). 

We would like to discuss our results for soft nonfactorizable contributions 
in comparison with the hard nonfactorizable effects calculated in the QCD- 
improved factorization approach. The best thing would be to calculate both soft, 
and hard contributions inside the same model. In principle, the light-cone sum 
rule approach presented here enables such a calculation, although the estimation 
of hard nonfactorizable contributions is technically very demanding, involving a 
calculation of two-loop diagrams. Therefore, we proceed with the QCD-improved 
factorization estimations for the hard nonfactorizable contributions. 
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We the inclusion of the hard nonfactorizable corrections, the ai parameter 
(35) looks as follows 



The estimations done the QCD-improved factorization [24] show hard-gluon 
exchange corrections to the naive factorization result of the order of ~ 25%, 
predicted by the LO calculation with the twist-2 kaon distribution amplitude. 
Unlikely large corrections are obtained by inclusion of the twist-3 kaon distri- 
bution amplitude. Anyhow, due to the obvious dominance of soft contributions 
to the twist-3 part of the hard corrections in the BBNS approach [11], it is 
very likely that some double counting of soft effects could appear if we naively 
compare the results. Therefore, taking only the twist-2 hard nonfactorizable cor- 
rections from [24] into account, recalculated at the ^ scale, our prediction (35) 
changes to 



The prediction still remains to be too small to explain the data. 

Nevertheless, there are several things which have to be stressed here in con- 
nection with the result. Soft nonfactorizable contributions are at least equally 
important as nonfactorizable contributions from the hard-gluon exchange, and 
can be even dominant. Soft nonfactorizable contributions are positive, and the 
same seems to be valid for hard corrections. While hard corrections have an imag- 
inary part, in the soft contributions the annihilation and the penguin topologies 
as potential sources for the appearance of an imaginary part were not discussed. 
A comparison between the result (37) and the experimental value \a^\ ~ 0.3 
for B — > J/'ijjK decay, with the recently deduced parameter from the color- 
suppressed B — > Z)M 7r decays, | «2 1 ~ 0.4 — 0.5 [10], provides clear evidence 
for the nonuniversality of the ai parameter in color-suppressed decays. 

5 Conclusions 

We have reviewed recent progress in the understanding of the underlaying dy- 
namics of exclusive nonleptonic decays, with the emphasis on the nonfactorizable 
corrections to the naive factorization approach. We have described in more detail 
QCD light-cone sum rule approach to the calculation of nonfactorizable contri- 
butions, and have shown results for B — > tttt [13] and B — > J/tpK [20] decays. 

The QCD-improved factorization method is reviewed in this volume by M. 
Neubert, [3]. 
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(36) 



a 2 = 0.16-0.19| 



(37) 
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